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ABSTRACT 

Differential astrometry measurements from the Palomar High-precision Astrometric Search for Ex- 
oplanet Systems have been combined with lower precision single-aperture measurements covering a 
much longer timespan (from eyepiece measurements, speckle interferometry, and adaptive optics) to 
determine improved visual orbits for 20 binary stars. In some cases, radial velocity observations exist 
to constrain the full three-dimensional orbit and determine component masses. The visual orbit of 
one of these binaries — a Com (HD 114378) — shows that the system is likely to have eclipses, despite 
its very long period of 26 years. The next eclipse is predicted to be within a week of 2015 January 24. 
Subject headings: astrometry - binaries: close - binaries: eclipsing - binaries:visual - tech- 
niques : interferometric 



1. INTRODUCTION 

Accurate measurements of fundamental properties of 
stars are necessary to develop models of stellar forma- 
tion, structure, and evolution. These are also key pa- 
rameters in understanding galactic mass-luminosity re- 
lationships. Combined visual-radial velocity (RV) orbits 
of binary stars provide measures of component masses 
and the distances to the systems (determining the intrin- 
sic stellar luminosities), while visual-only orbits coupled 
with parallax measurements can be used to measure the 
total mass of the system (see, fo r example IBoden et al.1 
[19991: IPourbaix fc Jorissenl |2000t Benedict et all 120011 ). 
Widely separated binaries are particularly important be- 
cause the stars are more likely to have evolved like single 
stars. The accuracies with which the stars' fundamental 
properties are known are improved by enhanced measure- 
ment precisions. Furthermore, establishing a set of high 
precision orbits will be useful as calibration standards 
for future p recision differential astrometry efforts be- 
ing pursued (lLazorenko et al.ll2009l iCameron et ali r2009: 
iHehniniak et all 120091 ) . 
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A technique has been developed to obtain high pre- 
cision (35 /xas) astrometry of close stellar pairs (sep- 
aration less than 1 arcsec; lLane fc Muterspaugh! 120041) 
using long-baseline infrared interf erometry at the Palo - 
mar Testbed Interferometer (PTI: iColavita et al.|[l999l) . 
This technique was applied to 51 binary systems as 
the Palomar High-precision Astrometric Search for Ex- 
oplanet Systems (PHASES) program during 2002-2008. 
PHASES science results included precision binary orbits 
and component masses, studies of the geometries and 
component physical properties in triple and quadruple 
star systems, and limits on the presence of giant planet 
companions to the binaries. 

Astrometric measurements were made at PTI, which 
was located on Palomar Mountain near San Diego, CA. 
It was developed by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology for NASA, as a testbed for 
interferometric techniques applicable to the Keck Inter- 
ferometer and other missions such as the Space Interfer- 
ometry Mission (SIM). It was operated in the J (1.2/im), 
H (1.6/xm), and K (2.2/^m) bands, and combined starlight 
from two out of three available 40-cm apertures. The 
apertures formed a triangle with one 110 and two 87 me- 
m % 1 ? J Bgglfir?e r s jn -feASES observations began in 2002 and 
continued through 2008 November when PTI ceased rou- 
tine operations. 

This paper is the second in a series analyzing the fi- 
nal results of the PHASES project as of its comple- 
tion in late 2008. The first paper describes the ob- 
serving method, sources of measurement uncertainties, 
limits of observing precisions, derives empirical scaling 
rules to account for noise sources beyond those predicted 
by the standard reduction algorithms, and presents the 
full catalog of astrometric m easurements from PHASES 
(Muterspaugh et al. 2010d). The current paper com- 
bines PHASES astrometry with astrometric measure- 
ments made by other methods as well as RV observa- 
tions (where available) to determine orbital solutions to 
the binaries' Keplerian motions, determining physical 
properties such as component masses and system dis- 
tance when possible. The third paper presents limits 
on the existence of substcllar tertiary companions orbit- 
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ing either the primary or secondary stars in those sys- 
tems that ar e found to be consistent w ith being sim- 
ple binaries (|Muterspaugh et all 1201 0b 1 ). The fourth 
paper presents three-component orbital solutions to a 
known triple star system (63 Gem = HD 58728) and 
a newly discovered t r iple sy stem (HR 2896 = HD 60318) 
(jMuterspaugh et all l2010al) . Finally the fifth paper 
presents candidate substellar comp anions to PHASES 
binari es as detected by astrometry (jMuterspaugh et all 
[20T0cT ). 

2. PHASES MEASUREMENTS 

PHASES differential astrometry measurements were 
obtained using the observing method and standard data 
analysis pipeline described in Paper I. The measure- 
ments themselves and associated measurement uncer- 
tainties are also tabulated in Paper I. Analysis of bi- 
nary orbits is limited to the 20 systems for which 10 or 
more PHASES observations were made, do not show con- 
vincing evidence of having tertiary companions, and are 
not 6 Equ, which has al ready been the sub ject of othe r 
PHASES investigations (jMuterspaugh et al.ll2005ll2008t ). 
Those having stellar companions to the primary and/or 
secondary stars either have been pr es ented in previous 
works (see iMuterspaugh eFa l. 2006b,a| : j Lane et al.ll2007t 
iMuterspaugh et al.ll2008t lLane et al.l " l2O10l ) or in Paper 
IV. Those for which the PHASES measurements indicate 
substellar companions may be present are the subject of 
Paper V. The number of PHASES measurements avail- 
able for each of the 20 systems being investigated are 
listed in Tableffl 

3. NON-PHASES MEASUREMENTS 

Measurements of binaries observed by PHASES made 
by previous astrometric techniques and cataloged in 
the Washington Do uble Star Catalog (WPS , and see 
references therein, IMason et all [200lL [2010) were as- 
signed weights accor ding to the formula described by 
lHartkopf et ail pool . Unit weight uncertainties in sep- 
aration and position angle were evaluated by the follow- 
ing iterative procedure. First, guess values for the unit 
uncertainties of 24 mas in separation and 1.8° in posi- 
tion angle were assigned to the measurements of a given 
binary; these values correspon ded to previous experienc e 
using this procedure on fi Ori (jMuterspaugh et al~l [2008'). 
Second, the measurements were fit to a Keplerian model 
and the orbital parameters were optimized to minimize 
the fit x 2 - Third, the weighted scatter of the residuals 
in separation and position angle were evaluated, and the 
guessed unit uncertainties updated to make the rms scat- 
ter in each equal to unity. Fourth, the second and third 
steps were iterated two more times, at which point the 
values converged. Fifth, the final unit uncertainties were 
multiplied by the square root of the reduced \ 2 (V^r) 
of the fit, and refit one more time using these slightly 
larger weights. Sixth, if no residuals deviated by more 
than 3cr, the process ended, otherwise, the single mea- 
surement with the largest separation or position angle 
residual (weighted by its uncertainty) was flagged as an 
outlier and removed from future fits. Seventh, the pro- 
cess was repeated at the first step. The resulting weights 
are listed in Table [T] and the measurements themselves 
are listed in Table [2j 



Table 1 

Number of PHASES and Non-PHASES Measurements and Unit 
weight Uncertainties for Non-PHASES Measurements 



HD Number 


N P 


Np,o 


Nnp Nnp,o 


°>,° 




5286 


18 





645 


34 


0.094 


1.82 


6811 


19 


1 


185 


18 


0.031 


2.86 


17904 


45 


1 


150 


20 


0.017 


3.95 


26690 


19 


1 


96 


16 


0.0088 


3.11 


44926 


23 





11 





0.023 


5.57 


76943 


16 





217 


9 


0.060 


4.17 


77327 


17 





125 


12 


0.018 


2.61 


81858 


12 


o 


580 


16 


0.090 


2.23 


114378 


23 


1 


586 


33 


0.068 


2.09 


129246 


16 


1 


682 


15 


0.092 


1.95 


137107 


50 


1 


948 


63 


0.075 


2.05 


137391 


22 


1 


12 





0.0055 


3.25 


137909 


73 





93 


13 


0.016 


2.96 


140159 


15 





205 


13 


0.021 


3.78 


140436 


37 


6 


409 


26 


0.078 


3.08 


155103 


10 





131 


23 


0.0093 


3.44 


187362 


10 





223 


23 


0.021 


2.90 


202444 


39 





286 


13 


0.123 


4.37 


207652 


50 





185 


11 


0.030 


4.06 


214850 


48 


5 


253 


21 


0.036 


3.18 



Note. — The number of PHASES and non-PHASES astromet- 
ric measurements used for orbit fitting with each of the binaries 
being studied are presented in Columns 2 and 4 respectively, along 
with the additional numbers of measurements rejected as outliers 
in Columns 3 and 5. Columns 6 and 7 list the la measurement un- 
certainties for unit weight measurements from non-PHASES obser- 
vations determined by iterating Keplerian fits to the measurements 
with removal of 3<r or greater outliers in either dimension. Columns 
6 and 7 are in units of arcseconds and degrees, respectively. The 
subscripts "P" and "NP" stand for PHASES and Non-PHASES 
type measurements, respectively, while the subscript "O" stands 
for outlier and the subscript "o" represents that the uncertainties 
are for unit- weight measurements. 

4. ORBITAL SOLUTIONS 

4.1. Binaries Without Radial Velocity Measurements 

Binaries for which only astrometric measurements are 
available were fit to a single Keplerian orbit using a down- 
hill \ 2 minimizing routine alternating between the stan- 
dard Campbell and Thiele-Innes parameter sets for im- 
proved speed of convergence. When only a visual or- 
bit is available, the sum of the component masses can 
be determined if a trigonometric parallax measurement 
is available. For each system evaluated, the trigono- 
metric parallax and its associated measurement uncer- 
tainty, as me a sured by the revised Hipparcos analysis of 
Ivan Leeuw cn (2007), were taken as input parameters to 
the fit in order to evaluate the mass sums. When avail- 
abl e, the revised par allax estimates based on the work 
of (ISoderhielml [l99l , hereafter S99) (HD 76943, HD 
114378, HD 137391, HD 140159, HD 140436, HD 155103, 
HD 202444, and HD 207652) were used instead of the re- 
vised Hipparcos values because the S99 analysis made 
more complete use of non- Hipparcos differential astrom- 
etry values to separate the binary orbit signal from that 
of parallax. 

For each of the 15 binaries, the PHASES and non- 
PHASES astrometric measurements, and Hipparcos- 
based parallax measurements were combined in a single 
orbital fit to evaluate the Campbell orbital parameters 
and system distance. Alternatively, the orbit semimajor 
axis was replaced with the total system mass; repeat- 
ing the fit with this substitution allowed a natural way 
of evaluating the uncertainties of the total mass. While 
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Table 2 

Non-PHASES Astrometric Measurements 



HD Number 


Date 


P 


e 


0> 


O0 


Weight 


Outlier 




(yr) 


(arcsec) 


(dog) 


(arcsec) 


(deg) 






5286 


1830.7400 


0.850 


305.00 


0.297 


5.77 


0.1 





5286 


1831.7900 


0.780 


308.70 


0.297 


5.77 


0.1 





5286 


1832.1400 


0.850 


307.80 


0.086 


1.67 


1.2 





5286 


1835.9200 


1.100 


315.70 


0.297 


5.77 


0.1 





214850 


2008.5377 


0.363 


77.10 


0.008 


0.72 


19.4 





214850 


2008.5460 


0.364 


77.00 


0.008 


0.72 


19.4 





214850 


2008.5460 


0.364 


77.10 


0.008 


0.72 


19.4 





214850 


2008.8875 


0.380 


77.90 


0.024 


2.10 


2.3 






Note. — Non-PHASES astrometric measurements from the WDS are listed with ltr measurements uncertainties, and weights. Column 
1 is the HD catalog number of the target star, Column 2 is the decimal year of the observation, Columns 3 and 4 are the separation in 
arcseconds and position angle in degrees, respectively, Columns 5 and 6 are the lcr uncertainties in the measured quantities from Columns 
3 and 4, Column 7 is the weight assigned to the measurement, and Column 8 is 1 if the measurement is a >3cr outlier and omitted from 
the fit, otherwise. (This table is available in its entirety in machine-readable and Virtual Observatory (VO) forms in the online journal. 
A portion is shown here for guidance regarding its form and content.) 



the system distance is entirely dependent on the parallax 
measurement, including this as an input value with its 
associated uncertainty was necessary to evaluate the re- 
sulting mass uncertainties properly at the conclusion of 
the fitting procedure. The best-fit parameters and their 
uncertainties are listed in Table 2J 

4.1.1. HD 5286 

HD 5286 (36 And, HR 258, HIP 4288, and WDS 
00550+2338) is a pair of subgiant stars with spectral 
types G6 and K6. The total system mass of 1.86 ± 
0.15 Mp, is in r easonable agreement with this classifica- 
tion (1CqxII2000D . 

4.1.2. HD 6811 

HD 6811 O And, 42 And, HR 335, HIP 5434, and WDS 
01095+4715) is a pair of B stars (B6IV and B9V). The 
total mass is 6.5 ± 2.8 M ; while not well constrained, 
this value is reasonable for these stars. 

4.1.3. HD 17904 

HD 17904 (20 Per, HR 855, HIP 13490, and 
WDS 02537+3 820) was rep orted to have a 1269 day 
subsystem by lAbt fe Levvf (119761). The companion 
was further inves t igated by I Scarf e fe Fekell (|1978l ) and 
iMorbev fc Grif fin (1987j) who found no evidence of such 
a companion, nor is such an object found by astrometry. 
The single-component velocities from those works were 
consistent with constant velocity, providing no constraint 
on the binary orbit, and were not used for the present 
orbit fitting. HD 17904 is a pair of mid F dwarfs; the 
total system mass of 3.88 ±0.58 Mq is a bit high for such 
stars, though the uncertainty in this is large enough to 
still be in reasonable agreement. 

4.1.4. HD 44926 

While the orbital parameters of HD 44926 (HIP 30569, 
WDS 06255+2327) are not well constrained due to few 
observations covering a short fraction of the orbit, this is 
the first time an orbital solution has been evaluated for 
which the fit converged. 

4.1.5. HD 76943 

iHeintd (|1981l ) published four single-component spec- 
troscopic velocities of HD 76943 (10 UMa-though it is 



now in the constellation Lynx ()Griffinl fl999t) . HR 3579, 
HIP 44248, and WDS 09006+4147); the first was taken a 
year before the others, which were all on the same night. 
More recently, Tennessee State U niversity's Automatic 
Spect roscopic Telescope (AST; lEaton fc Williamson! 
2007) collected 66 blended spectra of 10 UMa spanning 
1843 days (23 % of the binary orbit). Upon inspection, 
the single set of absorption lines appeared asymmetric 
due to the binarity, but remained blended by an amount 
substantial enough to cause systematic errors in attempts 
at component velocity determinations. Higher resolution 
observations at times of maximal velocity separation, or 
individual component spectroscopy behind an Adaptive 
Optics system may allow better modeling of each compo- 
nent's lines, in which case these measurements could be 
recovered, but at this time the system remains without 
useful velocity measurements. 

The spectral type is consistently listed as F5V, con- 
sistent with the total system mass of 2.44 ± 0.12 M Q 
(|Coxl I2000D and the S99 component mass values of 
1.37 ± 0.08 M Q and 1.04 ± 0.06 M ; 

4.1.6. HD 77327 

A few single component RV measurements of HD 77327 
(k UMa, 12 UMa, HR 3594, HIP 44471, and WDS 
090 36+470 9 ) hav e been published by lAbt et al.l (|1980f ) 
and IHeirit^ (|1981f) . However, these are small in number, 
span a short time period, and appear to be consistent 
with constant velocity. As a result, the system is treated 
as a visual binary only. 

HD 77327 is composed of a pair of early (~ AO) dwarf 
stars. The total mass of 6.30+0.98 M Q is not constrained 
very well, and while the value is a bit high for stars of this 
type, the large uncertainty indicates there is no cause for 
alarm. 

4.1.7. HD 114378 

CORAVEL produced 22 two component RV measure- 
ments of HD 114378 (a Com, 42 Com, HR 4968, HIP 
64241, and WDS 13100+1732) spanning nearly half 
(46%) of the binary orbit, making it bo th a visual binary 
and d ouble-lined spectroscopic binary (|Duauennov et all 
I1991D . However, despite the anticipated velocity ampli- 
tude of each star being ~ 8kms _1 -several times larger 
than the ~ lkms -1 typical measurement uncertainty- 
no significant velocity changes appear in the CORAVEL 
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velocities. As a result, the velocities did not agree well 
with the visual orbit (after orbit fitting, over half have 
residuals marking them as outliers of > 4er compared 
to the published measurement uncertainties), and pro- 
duced a system distance that was much too distant to 
agree with that of S99 (55.9 ±1.4 mas) based on Hippar- 
cos astrometry. Thus, these do not appear to be useful 
in orbit fitting. 

Alternatively, iTokovinin fc Smekhovl (|2002t ) measured 
nine single-component RVs spanning ~ 13% of the binary 
orbit that were consistent with the visual orbit, though 
these are clustered at two observing epochs and show- 
ing little variation, making spurious orbit fitting likely. 
Using these velocities and the parallax of 55.9 ± 1.4 mas 
from S99 yielded a result that is almost certainly spuri- 
ous; for example, the binary mass ratio was a very low 
value (0.02). Thus, only the visual orbit is presented at 
the current time, and additional investigations will need 
to concentrate on improving the velocity results. 

The total mass was found to be 2.45 ± 0.18 M©, near 
the value of 2.54 M© of S99, though this is not surprising 
given that the parallax values used were identical. The 
system is a pair of F5 dwarf stars, and this mass range 
is reasonable for that spectral type. 

Most interesting of all, it is likely that this long period 
(26 y ears) binary eclipses! This has been anticipated be- 
fore (jHartkopf et al.lll989HHoffleitlll996l) . but previously 
the orbit was not determined well enough to demonstrate 
that eclipses are likely. The visual orbit parameters from 
the present investigation predict a closest approach of 
0.32 mas around 2015 January 24 (to within about a 
week), when the star is observable during early morning 
hours, and lasting about 1.5 days. Mid-F dwarfs have 
physical sizes about 1.3 times larger than the Sun. At a 
distance of 17.9 pc, this corresponds to diameters of 0.7 
mas. Eclipses are likely even if one considers the la up- 
per limits of distance and orbital inclination, for which 
closest approach increases only to 0.38 mas and the ap- 
parent stellar diameters decrease by 2.5%, an amount 
much smaller than the uncertainty in the stars' physi- 
cal sizes; see Figure [1] Unfortunately, only one eclipse 
per orbit is likely for this pair; the eccentric orbit causes 
the alternative closest separation to be 0.9 mas, and a 
likely near-miss, a Com, the brightest star in its con- 
stellation, with an orbital period of 26 years (nearly as 
long as e Aur) is likely to join e Aur as representing long 
period eclipsing binaries, despite the low probability of 
such a configuration. Unlike e Aur, in which the chance 
of eclipse is enhanced by the presence of an extended 
disk around one star, the eclipses in a Com result purely 
from orbital geometry and the overlapping of the stars 
themselves. 

4.1.8. HD 129246 

HD 129246 « Boo, 30 Boo, HR 5477, HIP 71795, and 
WDS 1411+1344) has an extremely high eccentricity of 
0.9977 ±0.0034. The distance of closest approach is only 
0.3 AU. Projected on the sky, the distance of closest ap- 
proach is 1 mas; in this case, eclipses are unlikely, though 
the uncertainty on the inclination is large. The next time 
of closest approach is 2023 August, though at that time 
the star is not up at night. 

Though the system has a very long orbital period of 
124 years, nearly two orbits have been observed; the first 
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t (MJD) 

Figure 1. Time of the predicted eclipse of HD 114378 (a Com) 
in late 2015 January. The solid line shows the projected sky sep- 
aration of the stars in the binary as a function of time for the 
best-fit orbital solution. The thinner dashed lines show the sep- 
aration predicted by changing the value of the orbital inclination 
by ±1<t. Varying other elements of the orbital solution by their 
lcr uncertainties does not significantly change the angle of closest 
approach, but does change the predicted time of the eclipse. The 
horizontal line at 0.7 mas corresponds to the expected physical 
diameter of the stars. 

measurement was from 1796. Despite this time coverage, 
and as a result of the large eccentricity, the semimajor 
axis remains poorly determined. As a result, the total 
mass is found to be very large (122 Mq) but also very un- 
certain (±275 Mq) to the point of being unconstrained. 

4.1.9. HD 137391 

The AST attempted observations of HD 137391 
{H Boo, 51 Boo, HR 5733, HIP 75411, and WDS 
15245+3723) several times, but the absorption lines were 
always too broad and blended for RV measurements. The 
mass sum of 3.24 ± 0.23 M Q is smaller than the value of 
3.70 Mq from S99 despite using the same distance to 
the system. The current value of the semimajor axis is 
smaller than that from S99, explaining the discrepancy. 

4.1.10. HD 140159 

HD 140159 (t Ser, 21 Set, HR 5842, HIP 76852, and 
WDS 15416+1940) is a pair of early A dwarfs with mass 
3.73 ± 0.53 M Q , consistent with both the S99 value of 
3.79 Mq and their spe ctral types , though slightly on the 
low side for the latter (|Coxj|200Q[ ). 

4.1.11. HD 140436 

HD 140436 (7 CrB, 8 Crb, HR 5849, HIP 76952, and 
WDS 15427+2618) is a pair of early A stars. The total 
mass of 3.73 ± 0.37 Mq is slightly lower than the value 
of 4.23 Mq of S99, though still in reasonable agreement. 
Like HD 140159, this mass is a bit low for their spectral 
types. 



4.1.12. HD 155103 

With only 10 PHASES measurements being available, 
the orbit model of HD 155103 (c Her, HR 6377, HIP 
83838, and WDS 17080+3556) was not much altered by 
their addition. The total mass of 3.32 ± 0.39 Mq is iden- 
tical to that from S99. The spectral types are early (A 
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or F), though the spectral class has been reported as ei- 
ther dwarf or giant. The total mass appears to be more 
consistent with dwarf stars. 

4.1.13. HD 187362 

HD 187362 (C Sge, 8 Sge, HR 7546, HIP 97496, and 
WDS 19490+1909) is a pair of early A dwarfs, though 
the measured total mass of 2.23 ± 0.35 Mp, is a bi t lower 
than one would expect for such stars (jCoxl 12000). This 
may be due to a remaining error in the parallax: if the 
original Hipparcos parallax is used instead, the total mass 
is found to be 4.7 + 1.2 Mq, which seems more consistent 
(|Perrvman et al.lll997l ). 

4.1.14. HD 202444 

Though the binary period of HD 202444 (r Cyg, 65 
Cyg, HR 8130, HIP 104887, and WDS 21148+3803) is 
long and the velocity amplitudes are low, it is possible 
to use spectroscopy to determine individual velocities for 
both components. This is because one star is extremely 
broad lined while the other has very sharp spectral fea- 
tures. Though the spectral lines are always blended, this 
odd pairing makes a very distinctive line profile that 
can be fit a s the sum of two very different Gaussians. 
iFekel et all (|2003l ) published two two-component RVs of 
the system. The TSU AST recently observed r Cyg 
with high cadence to better constrain the mass ratio and 
search for short period companions to either star. No 
significant changes were observed in the velocities, and 
the span of observations was too short to contribute to 
orbital modeling, so at this time only the system's visual 
orbit was evaluated. 

There is some indication that r Cyg may have a sub- 
stellar companion orbiting one of the two stars (see Pa- 
per V). There are reasons to doubt the authenticity of 
this proposed companion, so the visual orbit obtained by 
modeling the system with only a single Keplerian model 
has been presented here in addition to the double Ke- 
plerian model presented in Paper V. If real, the com- 
panion has a long orbital period. This signal could be 
absorbed into that of the wider binary when only the 
shorter timespan PHASES data were analyzed to search 
for tertiary companions, so no compelling evidence for 
a companion was present when only PHASES measure- 
ments were analyzed. However, the continued large value 
of x 2 that results when the combined PHASES and non- 
PHASES astrometry set described in Section [3] was ana- 
lyzed prompted a second search for tertiary companions, 
this time using all the astrometric measurements. The 
longer timespan non-PHASES astrometry reduced the 
amount by which the binary orbit parameters could ab- 
sorb motion of an intermediate period companion (short 
compared to the binary motion, long compared to the 
timespan of PHASES measurements) and indicated the 
presence of a companion with mass corresponding to that 
of a giant planet. 

4.1.15. HD 207652 

HD 207652 (13 Peg, HR 8344, HIP 107788, V373 Peg, 
and WDS 21501+ 1717) is a variable flar e star. It has 
been suggested by iTamazian et al.l (|1999l ) that the sec- 
ondary may be a T Tauri type variable. However, the 
variability does not appear to have affected the astrom- 
etry measurements; see Figure [2j The total mass is 



2.65 ± 0.21 Mq, in good agreement with the value of 
2.67 Mq from S99. For the system's early F spectral 
type, this mass and the overall system luminosity are 
more consistent with dwarf stars, though the spectral 
classification is giant or subgiant. 

4.2. Binaries With Single- component Radial Velocity 
Measurements 

Spectroscopic RVs can improve orbit fitting and enable 
component masses to be measured. If only one star's 
spectral features can be observed (single-lined spectro- 
scopic binaries), the masses can be measured if the sys- 
tem's distance has been measured via trigonometric par- 
allax and if a sufficient part of the orbit has been studied 
spectroscopically. Such measurements are available in 
the literature for three PHASES binaries: 46 Tau (HD 
26690), uj Leo (HD 81858), and /3 CrB (HD 137909). Of 
these, a large number of measurements covering several 
orbital periods were available for ft CrB, while the oth- 
ers had only small amounts of data available spanning 
a fraction of the orbital period. For (3 CrB, including 
RV measurements aids constraining all aspects of the or- 
bit, whereas for 46 Tau and w Leo most orbital elements 
are unaffected, but the mass ratio (and thus individual 
masses) can be constrained (though at a limited level). 
A separate system velocity parameter V was used for 
each velocity data set of a given star to allow for instru- 
mental offsets (for example, there are six independent V Q 
parameters for /3 CrB, one for each paper from which ve- 
locity data were obtained). The results of these fits are 
presented in Table [5] 

4.2.1. HD 26690 

Five single-component velocity measurements covering 
488 days (18% of the binary orbit) of HD 26690 (46 Tau, 
HR 1309 HIP 19719, a nd WDS 04136+0743) were pub- 
lished bv iHeinta (|T981). An additional 26 unpublished 
measurements covering 256 days (10% of the binary or- 
bit) have been made by TSU's AST and are presented 
for the first time in Table |31 A combined fit between the 
astrometric and velocimetry measurements found typical 
uncertainties of 1.6 kms -1 and 0.65 kms -1 were needed 
for these data sets, respectively, to make the weighted 
rms of the velocimetry residuals equal to 1. The revised 
Hipparcos parallax by S99 of 27.41+0.93 mas was used as 
a weighted input in orbit fitting, as an observation with 
associated uncertainty. The combined fit solved for the 
system distance as a parameter, and both the direct par- 
allax measurement and orbital fit were used to constrain 
this in an optimal manner. The component masses are 
1.38 ± 0.26 M Q and 0.82 ± 0.21 Mq, a mass ratio larger 
than that from S99, in which masses of 1.15 ± O.14M 
and 1.10 ± 0.14M Q were found. 



4.2.2. HD 81858 

lAbt fc Levvl (|1976h published 20 spectroscopic mea- 
surements covering 1467 days (3.4% of the binary orbit) 
of HD 81858 (w Leo, 2 Leo, HR 3754, HIP 46454, and 
WDS 09285+0903). The combined astrometric and ve- 
locimetry orbital fit revealed that a measurement uncer- 
tainty of 1.0 kms -1 yielded a weighted rms of 1 for the 
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Figure 2. Fit residuals to HD 207652 (13 Peg) along the differential right ascension (left) and differential declination (right) axes for the 
PHASES measurements. The residuals have low scatter indicating that stellar variability did not affect the astrometry measurements. 



Table 3 

AST Velocities of 46 Tau 



Day (HMJD) 



RV (kms- 1 ) 



<7Ry(kms 1 ) 



55024.46106 


-4.16 


0.65 


55052.38144 


-4.70 


0.65 


55063.44779 


-5.13 


0.65 


55089.51270 


-5.60 


0.65 


55094.51579 


-5.02 


0.65 


55095.30008 


-4.34 


0.65 


55099.50654 


-4.48 


0.65 


55102.41085 


-4.81 


0.65 


55118.43687 


-3.66 


0.65 


55121.46221 


-5.93 


0.65 


55122.46229 


-5.94 


0.65 


55135.38774 


-5.35 


0.65 


55136.26277 


-4.51 


0.65 


55137.28781 


-4.90 


0.65 


55139.37828 


-4.93 


0.65 


55141.42839 


-4.72 


0.65 


55156.33458 


-5.30 


0.65 


55157.11192 


-5.57 


0.65 


55162.26850 


-5.83 


0.65 


55183.28776 


-3.83 


0.65 


55201.39165 


-7.01 


0.65 


55236.11569 


-5.27 


0.65 


55245.16062 


-5.58 


0.65 


55258.19066 


-5.43 


0.65 


55271.17012 


-5.70 


0.65 


55280.11226 


-5.55 


0.65 



Note. — Single-component RV measurements of 46 Tau (HD 
26690) from the AST. 

velocity residuals. The short span of the spectroscopic 
measurements compared to the orbital period resulted in 
only loose constraints being placed on the binary mass 
ratio and individual component masses. The S99 paral- 
lax of 27.5 ±1.1 mas was used as a weighted input in 
orbit fitting. The masses are only poorly constrained as 
1.92 ± 0.96 M and 0.28 ± 0.92 M ; RV measurements 
covering more of the binary orbit and/or double-line RV 
measurements are needed to improve these values. The 
total system mass could be better determined by im- 
proved parallax measurements, but additional RV mea- 
surements are necessary for improved evaluation of the 
mass ratio. 

4.2.3. HD 137909 

A considerable number of RV measurements of HD 
137909 ("Peculiar Rosette Stone", /3 CrB, 3 CrB, HR 
5747, HIP 75695, and WDS 15278+2906) have been pub- 



lished, spanning nearly a century. These include the fol- 
lowing: 



• 91 measurements by I Cannon! (|1912[ ) (a 
4.3 kms" 1 ), 



352 measurements by iNeubauerl (|1944f ) (a 
l^kms- 1 ), 



25 measurements bv lWolffl (|1978l ) (ct ~ 0.8 km s" 1 ) 
60 measurements bv lOetken fc Orwertl {1984) (cr " 



1.0 kms" 1 ), 

121 measurements by iKamper et al.l (|1990f ) (unit 
uncertainty a ~ 1.2 kms -1 ), and 



• 78 measur ements by INorth et al.l (|1998D ( uncer- 
tainties as given in that paper), 

where the measurement uncertainties in all but the last 
case were derived based on that which yields a weighted 
rms of 1 for the velocity residuals from a combined orbital 
fit to the astrometric and velocity data. The S99 parallax 
of 29.31 ± 0.82 mas was used as a weighted input in orbit 
fitting. 

The PHASES and non-PHASES astrometric measure- 
ments, velocity measurements, and S99 parallax mea- 
surement were combined in a single orbital fit to evaluate 
the Campbell orbital parameters as well as binary mass 
ratio and system distance. Alternatively, the orbit semi- 
major axis and mass ratio were replaced with either the 
total system mass and mass ratio or the masses of each 
component; repeating the fit with these substitutions al- 
lowed a natural way of evaluating the uncertainties of 
these degenerate quantities. While the system distance 
is entirely dependent on the parallax measurement, in- 
cluding this as an input value with its associated uncer- 
tainty was necessary to properly evaluate the resulting 
mass ratio and mass uncertainties at the conclusion of 
the fitting procedure. The masses are 1.71 ± 0.18 M Q 
and 1.330 ± 0.074 M^, values a little smaller than those 
from INorth et al.l (|1998l ). 

4.3. Double-lined Spectroscopic Binaries 

When the visual orbit of a binary as well as the RVs 
of both components are available, the component masses 
and distance to the system can be determined to very 
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high precision. The distance can be used to convert ap- 
parent brightnesses to absolute magnitudes, offering very 
strong constraints on the stars' physical properties. Two 
PHASES binaries have RV measurements available for 
both components and the resulting orbital solutions are 
presented here. 

One additional double-line spectroscopic binary ob- 
served with PHASES is 8 Equ (HD 202275). Nineteen 
new PHASES measurements of i5 Equ were made sinc e 
its orbit was last updated in iMuterspaugh et al.1 {2008), 
which included 49 measurements. However, none of the 
orbital elements, the component masses, the system dis- 
tance, nor any of their associated uncertainties, are sig- 
nificantly different when all 68 measurements were in- 
cluded as compared to that previous work, so it was not 
included in Table [6] because no updated information was 
needed. 

4.3.1. HD 137107 

At least one substellar companion has already been 
detected in the HD 137107 (77 CrB, 2 CrB, HR 5727, HIP 
75312, and WDS 15232+3017) system, though it has no 
impact on the PHASES study. iKirkpatrick et al.1 (|2001l ) 
imaged a faint L8 dwarf of mass 0.060 ± 0.15 M at 194" 
southeast of the A-B pair. They concluded that the pair 
is physically associated, with a separation of ~ 3600 AU. 
It is a circumbinary companion to the A-B binary. 

CORAVEL produced 31 two-component velocity mea- 
surements of HP 1371 07 A and B, reported first b y 
IDuauennov et all (|1991f) and revised bv lPourbaixl(|2000D . 
The TSU AST observed this system 33 times, but unfor- 
tunately these observations took place at a time when the 
lines were blended, and thus are not used for orbit fitting. 
The CORA VEL measurements have been combined with 
the astrometry to develop a full three-dimensional orbit. 
The S99 parallax of 53.5 ± 0.9 mas was also used as an 
additional measurement with associated uncertainty. 

The masses are 1.243±0.054M Q and 1.100±0.039 M Q . 
As reported in Paper I, the differential magnitude of rj 
CrB was measured in the K band using Keck Adaptive 
Optics. The differential magnitude is AKn = 0.185 ± 
0.001. The apparent K magnitude from iSkrutskie et al.1 
( 2006) is 3.714 ±0.216. Combined with the orbit-derived 
distance of 18.50 ± 0.22 pc (which is consistent with 
the trigonometric parallax) , these correspond to absolute 
magnitudes of K x = 3.04 ± 0.22 and K 2 = 3.23 ± 0.22. 

4.3.2. HD 214850 

iBatten et al.l (|1985| ) published 10 single-component ve- 
locities and 12 two-component velocity measurements 
of HD 214850 (HR 8631, HIP 111974, and WDS 
22409+1433) spanning 1614 days (~ 21% of the binary 
period) and 271 days (~ 3.6% of the binary period), 
respectively. Four more measurements are omitted, as 
done by the original authors, due to partial blending of 
the spectral lines. Unit uncertainties of 0.19kms _1 for 
component A, and 0.72 kins -1 for component B were as- 
signed to these measurements, resulting in the weighted 
rms scatters of each being unity after orbit fitting. The 
S99 parallax of 29.5 ± 0.8 mas is consistent with the or- 
bital solution and has been included as an input con- 
straint for orbit fitting and is reflected in the orbital pa- 
rameters and their uncertainties in Table |6l Component 
masses are measured at the 7% level. 



5. CONCLUSIONS 

The high-precision nature of the PHASES measure- 
ments and the long timespan of the previous astrometric 
measurements that aid in closing the orbits determine 
the visual orbits of many long period binaries with high 
accuracies. These orbital solutions will benefit future 
astrometric programs by acting as calibration sources. 
Future astrometric studies of these systems will aid the 
orbits presented by extending the time baseline of the 
precision observations. U nfortunately, it is unlikely the 
European Gaia mission (|Lindegren fc Perrvmanl Il996f) 
will contribute to the study of these specific systems be- 
cause its detectors saturate on bright stars, though it 
is possible the proposed J-MAPS mission could do so 
(jGaume fc Hen nessv 20091). In most cases, the stellar as- 
trophysics applications are limited by the RV data avail- 
able. In some cases, the lack of RV could be addressed 
simply by pursuing these measurements over the next 
decade or so. 
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Table 4 

Visual Orbit Parameters 



HD Number 


Period (d) 


T (HMJD) 


Scmimajor Axis (arcscc) 


Eccentricity 


Inclination (deg) 


lu (deg) 


U (deg) 


+ M 2 (M ) x 2 




<7p 


°"T 




0" e 


<7 i 




& il 


O Mi + M 2 


dof 


5286 


61183 


35543 


0.9837 


0.30603 


44.57 


358.62 


173.66 


1.86 


1574.1 




(69) 


(21) 


(0.0011) 


(0.00078) 


(0.11) 


(0.21) 


(0.13) 


(0.15) 


(1319) 


6811 


202459 


17740 


0.573 


0.385 


142.2 


112.6 


337.2 


6.5 


464.1 




(24530) 


(1837) 


(0.051) 


(0.043) 


(2.8) 


(9.1) 


(3.2) 


(2.8) 


(401) 


17904 


11553.9 


50255 


0.2224 


0.7560 


120.48 


265.54 


26.62 


3.88 


360.7 




(8.7) 


(12) 


(0.0011) 


(0.0023) 


(0.20) 


(0.11) 


(0.24) 


(0.58) 


(383) 


44926 


339008 


53508 


0.5 


0.7 


80.7 


14 


342.18 


2.53 


43.1 




(3132095) 


(2553) 


(3.0) 


(1.6) 


(4.9) 


(23) 


(0.70) 


(1.3) 


(67) 


76943 


7691.0 


49262.6 


0.64566 


0.15075 


131.366 


32.30 


203.74 


2.44 


476.9 




(1.8) 


(9.1) 


(0.00065) 


(0.00084) 


(0.099) 


(0.44) 


(0.10) 


(0.12) 


(459) 


77327 


13007.2 


50404 


0.18194 


0.5584 


109.410 


355.63 


105.641 


6.30 


297.0 




(9.7) 


(12) 


(0.00025) 


(0.0015) 


(0.066) 


(0.36) 


(0.080) 


(0.98) 


337 


114378 


9485.68 


47651.8 


0.66132 


0.4957 


90.054 


101.689 


12.221 


2.45 


1162.3 




(0.97) 


(2.6) 


(0.00061) 


(0.0010) 


(0.010) 


(0.059) 


(0.015) 


(0.18) 


(1211) 


129246 


45460 


60183 


2.3 


0.9977 


102.3 


262.9 


8.2 


122 


1452.2 




(62) 


(57) 


(1.7) 


(0.0034) 


(9.2) 


(5.9) 


(2.6) 


(275) 


(1389) 


137391 


1368.02 


53855.92 


0.098837 


0.27194 


130.016 


44.204 


130.040 


3.24 


52.1 




(0.24) 


(0.23) 


(0.000034) 


(0.00028) 


(0.029) 


(0.068) 


(0.040) 


(0.23) 


(61) 


140159 


8015.0 


54180 


0.21033 


0.0941 


83.608 


80.5 


69.684 


3.73 


453.1 




(6.3) 


(33) 


(0.00047) 


(0.0028) 


(0.043) 


(1.8) 


(0.033) 


(0.53) 


(433) 


140436 


33701 


59947 


0.7315 


0.4779 


94.263 


283.40 


292 043 


3.73 


1017.4 




(61) 


(63) 


(0.0034) 


(0.0067) 


(0.028) 


(0.36) 


(0.036) 


(0.37) 


(885) 


155103 


2971.54 


42567.9 


0.10863 


0.5338 


121.221 


53.81 


309.30 


3.32 


274.3 




(0.84) 


(3.4) 


(0.00020) 


(0.0015) 


(0.077) 


(0.16) 


(0.11) 


(0.39) 


(275) 


187362 


8487.9 


44199.6 


0.13605 


0.7948 


132.33 


355.3 


340.97 


2.23 


467.3 




(4.9) 


(5.9) 


(0.00044) 


(0.0019) 


(0.41) 


(1.0) 


(0.68) 


(0.35) 


(459) 


202444 


18122.6 


47527 


0.9160 


0.2413 


134.11 


298.28 


339.72 


2.63 






(8.8) 


(18) 


(0.0013) 


(0.0013) 


(0.14) 


(0.20) 


(0.14) 


(0.12) 


(643) 


207652 


9622.8 


47848.3 


0.36245 


0.2300 


70.379 


250.23 


230.698 


2.65 


431.8 




(4.1) 


(6.2) 


(0.00026) 


(0.0011) 


(0.029) 


(0.14) 


(0.028) 


(0.21) 


(463) 


IN UIK, 


The model parameters and fit uncertainties derived from a simultaneous fit to PHASES and non- 


PHASES measurements. 














Table 5 


















Single-line Spectroscopic Binary Orbit Parameters 






HD Number 


Period (d) 


To (HMJD) 


Semimajor Axis (arcscc) 


Eccentricity 


Inclination (deg) 


u (deg) 


n (deg) 


M 2 /Mi Mi 


+ M 2 (M Q 




(Tp 


(TT a 










en 


cr M 1 /M 2 


CM1+M2 


26690 


2633.34 


53415.37 


0.13302 


0.3301 


66.942 


306.658 


144.296 


0.60 


2.20 




(0.50) 


(0.39) 


(0.00015) 


(0.0010) 


(0.039) 


(0.071) 


(0.040) 


(0.23) 


(0.23) 


81858 


43089 


36649 


0.8599 


0.5619 


65.347 


302.54 


325.982 


0.14 


on 
z.zu 




(26) 


(11) 


(0.0022) 


(0.0014) 


(0.093) 


(0.14) 


(0.074) 


(0.55) 


(0.27) 


137909 


3848.54 


44412.8 


0.204008 


0.53971 


111.452 


180.21 


148.041 


0.779 


3.04 




(0.52) 


(1.4) 


(0.000034) 


(0.00021) 


(0.014) 


(0.13) 


(0.030) 


(0.039) 


(0.25) 


j\T te. 


The model parameters and fit uncertainties derived from a simultaneous fit to PHASES and non- 


PHASES measurements. 














Table 6 


















Double-line Spectroscopic Binary Orbit Parameters 




HD Number 


Period (d) 


To (HMJD) 


Semimajor Axis (arcscc) 


Eccentricity 


Inclination (deg) 


oj (deg) 


Q (deg) 


M2/M1 Mi + M 2 (Mq 




(Tp 


ff T„ 


c a 










<?M 1 /M 2 


0Mi+M 2 


137107 


15204.9 


42612.9 


0.86226 


0.27907 


58.084 


39.885 


202.827 


0.885 


2.343 




(1.4) 


(3.4) 


(0.00033) 


(0.00026) 


(0.026) 


(0.064) 


(0.024) 


(0.031) 


(0.084) 


214850 


7607.7 


45531.7 


0.287980 


0.73499 


139.861 


22.31 


251.540 


1.089 


2.246 




(1.1) 


(1.2) 


(0.000049) 


(0.00014) 


(0.032) 


(0.12) 


(0.076) 


(0.080) 


(0.067) 



Note. — The model parameters and fit uncertainties derived from a simultaneous fit to PHASES and non-PHASES astromctry and 
two component RV measurements. 



